Giant topological insulator gap in graphene with 5d adatoms 
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Two-dimensional topological insulators (2D TIs) have been proposed as platforms for many intriguing appli- 
cations, ranging from spintronics to topological quantum information processing. Realizing this potential will 
likely be facilitated by the discovery of new, easily manufactured materials in this class. With this goal in mind 
we introduce a new framework for engineering a 2D TI by hybridizing graphene with impurity bands arising 
from heavy adatoms possessing partially filled d-shells, in particular osmium and iridium. First principles cal- 
culations predict that the gaps generated by this means exceed 0.2 eV over a broad range of adatom coverage; 
moreover, tuning of the Fermi level is not required to enter the TI state. The mechanism at work is expected to 
be rather general and may open the door to designing new TI phases in many materials. 

PACS numbers: 



Topological insulators comprise a class of strongly spin- 
orbit-coupled, non-magnetic materials that are electrically in- 
ert in the bulk yet possess protected metallic states at their 
boundary [QJO. These systems are promising sources for a 
host of exotic phenomena — including Majorana fermions 0- 
7], charge fractionalization [8], and novel magneto-electric 
effects I9UT21 — and may also find use for quantum comput- 
ing (2) and spintronics devices 1 13 ). In some respects two- 
dimensional (2D) topological insulators are ideally suited for 
such applications; for example, bulk carriers that often plague 
their three-dimensional counterparts can be vacated simply by 
gating. Experimental progress on 2D topological insulators 
has steadily advanced recently due largely to pioneering work 
on HgTe lfT4lfT7l (see also Ref. 1 18]). Nevertheless, to real- 
ize their full potential systems more amenable to experimental 
investigations are highly desirable. In this regard the ability 
to design new 2D topological insulators from conventional, 
widely available materials would constitute a major step for- 
ward, and many proposals of this spirit now exist lfT9l426l . 

Following this strategy, here we introduce a new mech- 
anism for engineering a topological insulator state in 
graphene — arguably now the most broadly accessible 2D 
electron system. Historically, graphene was the first material 
predicted to realize a topological insulator in seminal work by 
Kane and Mele [ 1 ], though unfortunately the gap is unobserv- 
ably small due to carbon's exceedingly weak spin-orbit cou- 
pling I27U3T1 . Reference l23l revived graphene as a viable 
topological insulator candidate by predicting that dilute con- 
centrations of heavy In or TI adatoms dramatically enhance 
the gap to detectable values of order 0.01 eV. Essentially, these 
adatoms mediate enhanced spin-orbit interactions of the type 
present in the Kane-Mele model Q for pure graphene. 

Our approach here also relies on hybridizing graphene with 
dilute heavy adatoms, though the underlying physics is en- 
tirely different and can not be understood in terms of an effec- 
tive graphene-only model. Rather, we will show using density 
functional theory that certain adatoms — specifically Os, Ir, 
Cu-Os dimers, and Cu-Ir dimers — form spin-orbit- split impu- 
rity bands that hybridize with graphene 's Dirac states in such a 
way that a highly robust topological insulator regime emerges. 



In fact here it is more appropriate to view the adatoms as the 
dominant low-energy degrees of freedom, with their coupling 
effectively mediated by graphene; from this perspective this 
mechanism represents the inverse of that invoked in Ref. l23l . 

Numerous practical advantages arise in this scheme. The 
topological insulator gaps are extremely large — typically ex- 
ceeding 0.2 eV — and take on nearly the full atomic spin-orbit 
splitting for the adatoms. Such values reflect more than an 
order-of-magnitude enhancement compared to the gaps in- 
duced by In or TI, and are competitive with the largest gap 
predicted for any topological insulator. Somewhat counterin- 
tuitively, these gaps are remarkably insensitive to the adatom 
concentration, taking on comparable values at least over cov- 
erages ranging from ~ 2% — 6%. In the case of Os adatoms 
and Cu-Ir dimers, the Fermi level also naturally resides within 
the topological insulator band gap. This eliminates perhaps 
the most serious challenge with In and TI, both of which sub- 
stantially electron-dope graphene even at quite low coverages. 
These features suggest that the observation of a topological in- 
sulator state in graphene may be within reach. 

We first elucidate the mechanism uncovered here using a 
tight-binding model that exposes the physics in a very trans- 
parent manner. Consider bd adatoms residing at positions 
R located at 'hollow' (H) sites in graphene as in Fig. [TJah 
For simplicity we retain only the d xz and d yz adatom states 
since these comprise the most important orbitals in our first- 
principles calculations. (Recall that d xz / yz orbitals arise from 
L z orbital angular momentum m = ±1 states.) We then 
model the composite system by a Hamiltonian H = H g + 
H a + H c [23]. The first term allows nearest-neighbor hop- 
ping for graphene: 



H 9 = ~ t Yl ( C t a Cr>a + H.C.), (1) 

a=t,l <rr') 



where c\ a adds an electron with spin a to honeycomb site r. 
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FIG. 1: (a) 4 x 4 supercell employed to simulate periodic H-site 
adatoms (cyan) at 6.25% coverage, (b)-(d) Corresponding tight- 
binding band structures calculated with graphene hopping strength 
t = 2.7 eV and adatom on-site energy e = —0.5 eV. The adatom- 
graphene hopping t c and adatom spin-orbit coupling A so are given 
by (b) t c = Asa = 0, (c) t c = 1.5 eV, A so = 0, and (d) t c = 1.5 eV, 
A so — 0.2 eV. When the Fermi level sits at the dashed green line, 
A so 7^ generates a topological insulator gap given approximately 
by the atomic adatom spin-orbit splitting. 



The second encodes couplings for the adatoms, 

Ha — e /mRa/mRa 
R L a=t,|m=±l 

+ ^2 A so(fi na S Z a(3 fltl(3 - /iiRa^Z-lR/S) ( 2 ) 

Here /^ Rq; fills the adatom <i-orbital at position R with mag- 
netic quantum number m = ±1 and spin a, e sets the or- 
bital energies relative to graphene's Dirac points, A so rep- 
resents spin-orbit coupling, and s z is a Pauli matrix. Fi- 
nally, H c hybridizes the adatoms with graphene. To ex- 
press this term it is convenient to define vectors ej that point 
from an adatom to the six surrounding carbon sites [see Fig. 
pTa)]. One can then construct linear combinations C m R = 



and write (23 



\m(j-l) 



CR +ej that carry angular momentum m 



H c = ~tcJ2J2 E ( iC! Lnaf m n a + H.c). (3) 

R a=^,lm=±l 

Let us now specialize to a periodic adatom arrangement 
characterized by the 4x4 supercell shown in Fig.[lja), with 
one adatom per cell (this corresponds to 6.25% coverage). 
Figure [TJb) illustrates the band structure with t = 2.7 eV, 
e = —0.5 eV, and t c = A so = 0. In this limit the adatoms 
produce a four-fold degenerate flat band, reflecting spin and 



orbital degeneracy. For the following discussion the precise 
location of these adatom states is unimportant, provided they 
intersect the carbon bands within ~ 1 eV of the Dirac points. 
Incorporating tunneling between the adatoms and graphene 
causes the flat bands to disperse, as shown in Fig. [TJc) for 
t c = 1.5 eV. Suppose now that the Fermi level resides at the 
dashed line in Fig.[TJc). Although the spectrum here exhibits 
a sizable energy gap near the K point, the system remains 
metallic due to band touchings at the zone center. The gapless 
excitations at zero momentum exhibit the following two cru- 
cial properties: (i) they arise from weakly perturbed adatom 
orbitals since at the zone center the nearest carbon bands re- 
side well over 1 eV away in energy, and (ii) they are pro- 
tected by time-reversal, spatial rotation, and SU(2) spin sym- 
metries that coexist when A so = 0. Breaking the last of these 
symmetries by turning on spin-orbit coupling thus produces 
a bulk gap given nearly by the atomic spin-orbit splitting for 
the adatoms, despite their dilute coverage. This key point is 
demonstrated in Fig. [TJd) for A so = 0.2 eV, which yields a 
0.32 eV gap that constitutes 80% of a single adatom's spin- 
orbit splitting. 

The gap opening indeed drives the system into a topological 
insulator phase. Since our Hamiltonian is inversion symmet- 
ric this can be readily verified by computing Fu and Kane's 
formula for the Z 2 invariant in Ref. [32). For additional ev- 
idence the solid curve in Fig. |2ja) plots the density of states 
(DOS) near zero energy for the same periodic adatom cover- 
age on a graphene strip with armchair edges along x and peri- 
odic boundary conditions along y. (Our strip consists of 128 
zig-zag 'rows' of carbon sites, with 80 sites per row.) Edge 
states characteristic of the topological phase produce a finite 
DOS inside of the bulk gap, and are clearly resolvable in the 
system size simulated over an energy window of 0.31 eV. As 
an example, an edge state with mid-gap energy E = 0.004 
eV appears in Fig. [2jb), where circles indicate adatom loca- 
tions while the shading represents the probability amplitude 
extracted from the wavef unction. 

Remarkably, the formation of a topological insulator by no 
means requires the periodic arrangements considered so far. 
In fact similar physics arises even for completely randomly 
distributed H-site adatoms. The dashed curve in Fig. [2jb) il- 
lustrates the DOS for the random case, again at 6.25% cov- 
erage. Even in our finite system one can easily resolve edge 
modes within a 0.21 eV energy range that is comparable to 
the bulk gap for the periodic case; see, e.g., the mid-gap state 
with energy E = 0.003 eV plotted in Fig.^c). 

Next we demonstrate using density functional theory (DFT) 
that the mechanism described above can be realized in 
graphene with 5d adatoms, notably Os and Ir. All DFT calcu- 
lations were carried out with the Vienna ab-initio simulation 
package (VASP) (33] [34) at the level of the local density ap- 
proximation (LDA) (35), including spin-orbit coupling unless 
specified otherwise. Positions of all atoms were fully relaxed 
using the conjugated gradient method for energy minimization 
until the calculated force on each atom became smaller than 
0.01 eV/A. Most results were obtained using the supercell in 



FIG. 2: (a) Density of states for periodic (solid curve) and random 
(dashed curve) adatoms at 6.25% coverage on a graphene strip with 
armchair edges along x and periodic boundary conditions along y. 
Parameters are the same as for Fig.[T|c). The finite density of states 
within the bulk gap reflects edge states, which remarkably survive 
even for randomly distributed adatoms. Examples of edge states for 
the periodic and random cases respectively appear in (b) and (c). 



Fig.[TJa), with one adatom per cell. For additional details see 
l43l . Below we report our results for Os and then turn to Ir. 

Since the thermal stability of adsorption structures is rele- 
vant for both experiments and applications, candidate adatoms 
should ideally exhibit large H-site binding energies [defined as 
Eh = i^(graphene) +i?(adatom) — £^(adatom/graphene)] 
and high diffusion energy barriers [defined as = 
^(Transition state) — ^(Ground state)]. Osmium satis- 
fies both criteria. The binding energy for Os at the H site 
in graphene is 2.42 eV — much larger than the 'top' (directly 
above a C) and 'bridge' (above the midpoint of a C-C bond) 
configurations for which is 1.70 and 1.59 eV, respectively. 
Moreover, the calculated diffusion barrier for an Os adatom to 
diffuse from an H site through the top site is found to equal 
the difference between the binding energies at these positions, 
0.72 eV. The barrier for diffusion through the bridge site is 
similarly given by 0.83 eV. Therefore Os adatoms should be 
stable over H sites even at room temperature. By contrast most 
3d transition metals have ~ 1 eV [36 ] and are much more 
mobile (36H37); for example, the diffusion barrier is only 0.40 
eV for Co on graphene (38). 

Figure |3ja) displays the DFT band structure for periodic H- 
site Os adatoms on graphene at 6.25% coverage using the 4 x 4 
supercell in Fig. [TJa)- Each Os adatom forms a charge state 
of +0.55e (based on the Bader charge division scheme [39 ]), 
indicating that the Os-graphene bonds mix covalent and ionic 
features. Clearly these bonds dramatically modify the char- 
acteristic Dirac bands at the K point of pure graphene simi- 
lar to Fig. [T] Most importantly, Os induces a large band gap 
A so = 0.27 eV, right at the Fermi level given by the green 
dashed line in Fig.[3ja). As in our tight-binding model the gap 
here results solely from spin-orbit coupling. (Without spin- 
orbit interactions a gapless spectrum arises; see the Supple- 
mental Material.) More precisely, the partial density of states 
(PDOS) for the Os bd orbitals displayed in Fig.[3jb) indicates 
that the gap arises from the hybridization between graphene' s 
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FIG. 3: (a) First principles band structure for Os on graphene at 
6.25% coverage. The green dashed line indicates the Fermi level, (b) 
Corresponding partial density of states (PDOS) for the Os 5d levels. 
The large gap A so visible in (a) arises from hybridization between 
graphene and the spin-orbit- split d xz / yz orbitals, as in our tight- 
binding model, (c) Coverage dependence of the gap for graphene 
with Os adatoms (circles) and Cu-Os dimers (triangles). 



7r states and the spin-orbit- split d xz and d yz adatom orbitals, 
also as in our tight-binding model. The PDOS for the d z 2, 
d x 2_ y 2, and d xy orbitals, by contrast, is concentrated at much 
lower energies. Thus the gap-opening mechanism introduced 
earlier indeed appears in the realistic Os/graphene system. 

The Os-induced gap depends exceptionally weakly on cov- 
erage. To illustrate this important point we performed sim- 
ulations of graphene with one Os adatom in 5 x 5, 7 x 7, 
and 10 x 10 supercells (corresponding to coverages of 4%, 
2.04%, and 1%). Circles in Fig.^c) show the DFT-predicted 
gaps, which remain close to 0.2 eV even at 1% coverage. This 
striking feature is actually rather natural since the local atomic 
spin-orbit splitting for the Os d xz and d yz orbitals essentially 
sets Aso- 

Strictly speaking, a true topological insulator phase does 
not arise in the DFT simulations described above since Os 
forms small spin and orbital magnetic moments of 0.45 (1b 
and 0.05 /ig, respectively. This produces visible splittings 
of the bands at the Y and M points corresponding to time- 
reversal-invariant momenta; see Fig. [3ja). One should keep 
in mind, however, that DFT can sometimes overestimate mo- 
ment formation. Nonetheless, even if a moment M s indeed 
appears in an experiment, there are practical means by which 
this can be quenched to zero to reveal a bona fide topologi- 
cal phase [44 j. One effective approach is to apply an exter- 
nal electric field e. Figure |4ja) illustrates that M s of Os on 
graphene depends sensitively on e. In particular e < trans- 
fers additional charge from Os to graphene and kills the mo- 
ment for e < —0.3 V/A. The electric fields required to restore 
time-reversal symmetry only weakly affect the band structure. 
See, for example, Fig.|4jb) corresponding to 6.25% Os cov- 
erage with e = —0.5 V/A, where a time-reversal-invariant 
topological insulator appears with a gap A so = 0.26 eV. 
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Co-doping provides another means to quench the Os mag- 
netic moment. To preserve the main features of the band struc- 
ture while attracting charge away from Os (as accomplished 
by a negative e), co-adsorbates should interact weakly with 
graphene and exhibit larger electronegativity than Os. Fol- 
lowing this guidance, we considered Cu, Ag, and Au in sev- 
eral configurations as described in the Supplemental Material. 
Whereas Os repels Ag and Au adatoms, Cu prefers to climb 
over Os to form a vertical Cu-Os dimer over the H site. The 
binding energy = E (graphene) + E{Cn) + E(Os) — 
E(Cu — Os/graphene) for these dimers is 5.96 eV, higher 
by 2.50 eV compared to that of well- separated Cu and Os 
adatoms. Additionally, Cu more strongly anchors Os to the 
H site since the binding energy for the vertical dimers over 
the top (bridge) site is weaker by 1.27 (1.42) eV. This is not 
the only relevant criterion, however, since in practice isolated 
Cu and Os adatoms must be able to dimerize without over- 
coming substantial energy barriers. We explored this issue by 
computing the total energies along the diffusion path depicted 
in Fig.|4jc), where a Cu adatom beginning at position D ends 
up above an Os at position A (45). Figure |4jd) illustrates the 
change in energy AE relative to the dimer state for various 
positions along this trajectory. The energy barrier for a Cu 
adatom to diffuse from location D to B is only ~ 0.08 eV; 
once at position B the Cu strongly attracts to the top of Os 
and forms the vertical dimer l46l . This suggests that dimer 
formation ought to proceed quite efficiently. 

Because of the hybridization and charge transfer between 
the Cu and Os atoms — the Bader charges of Cu and Os are re- 
spectively — 0.21e and +0.67e — DFT predicts that graphene 
with Cu-Os dimers is nonmagnetic. The spectrum for (Cu- 
Os)/graphene at 6.25% coverage again supports a large topo- 
logical insulator gap A so = 0.21 eV as evident in the band 
structure of Fig. |4j e )- Moreover, the triangles in Fig. |3jc) 
show that this gap exhibits similarly weak coverage depen- 
dence as for Os/graphene. The drawback here, however, is 
that the Fermi level [green line in Fig. |4j e )] now resides in 
the valence band. Returning the Fermi level to the insulat- 
ing regime should be possible with conventional gating tech- 
niques, provided one works at low coverage. 

Alternatively, the hole introduced by each Cu-Os dimer can 
be compensated by replacing Os with Ir, which has one ad- 
ditional electron. Our calculations show that vertical Cu-Ir 
dimers also strongly bind to the H-site in graphene without 
forming a magnetic moment. Hybridization between Cu-Ir 
dimers and graphene produces nearly the same band structure 
as for (Cu-Os)/graphene, but with the Fermi level lying in the 
band gap. See the band structure for 6.25% Cu-Ir coverage in 
Fig.[4]X), where the gap is A$o = 0.25 eV. Additional results 
for Ir/graphene — which behaves similarly to Os/graphene — 
can be found in the Supplemental Material. 

In summary, we have introduced a mechanism by which 
graphene covered with heavy adatoms realizes a topological 
insulator protected by a giant gap comparable to atomic spin- 
orbit energies, even at exceptionally dilute coverages. Us- 
ing DFT we predicted that Os, Ir, Cu-Os dimers, and Cu-Ir 
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FIG. 4: (a) Magnetic moment M s of Os/graphene versus external 
electric field e applied perpendicular to the graphene sheet (see in- 
set for the direction of positive and negative e). The moment dis- 
appears for e < —0.3 V/A, allowing a true topological insulator 
phase to appear, (b) Band structure of Os/graphene with e = —0.5 
V/A. Comparing to Fig. [3ja) one sees that the fields needed to re- 
store time-reversal symmetry modify the band structure very little; 
in particular a large gap remains at the Fermi level (green dashed 
line), (c) Possible diffusion path of a Cu atom beginning from posi- 
tion D and ending above an Os atom at position A. (d) Energy profile 
along this diffusion trajectory. The small diffusion barrier evident 
in (d) indicates that Cu-Os dimers should form quite efficiently, (e) 
Band structure for Cu-Os dimers on graphene at 6.25% coverage. 
Time-reversal symmetry is preserved here even at s — 0, though the 
Fermi level now resides in the valence band, (f) Band structure for 
Cu-Ir dimers on graphene at 6.25% coverage. This system preserves 
time-reversal symmetry, eliminates the shift in Fermi level, and also 
supports a large topological insulator gap. 



dimers all give rise to this mechanism and produce gaps ex- 
ceeding 0.2 eV at coverages as low as 2%. Although our 
DFT calculations of necessity invoked periodic adatom con- 
figurations, our tight-binding simulations indicate that read- 
ily observable bulk (mobility) gaps should survive also in the 
random case relevant for experiments. These findings are ex- 
pected to greatly facilitate the realization of a topological insu- 
lator phase in graphene-based systems. We suspect, however, 
that the mechanism exposed here has much broader applica- 
tions since (contrary to Ref. [23]) the physics has nothing to 
do with the Kane-Mele model specific to graphene. Hybridiz- 
ing trivial metals or insulators with heavy-element impurity 
bands may therefore provide a generic method for designing 
new topological phases, which would be interesting to inves- 
tigate in future work. 

The authors gratefully acknowledge A. Damascelli, J. 
Eisenstein, J. Folk, E. Henriksen, and C. Zeng for helpful dis- 
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FIG. 6: Different configurations for Cu-Os co-adsorption. The dark 
cyan, coral, and gray spheres represent Os, Cu, and C atoms, respec- 
tively. 



Os/graphene and (Cu-Os)/graphene the gaps indeed originate 
solely from spin-orbit interactions. 



FIG. 5: Electronic properties of Os/graphene and (Cu-Os)/graphene 
at 6.25% coverage without spin-orbit coupling, (a) Partial density 
of states (PDOS) for the Os 5d orbitals in Os/graphene. Positive 
and negative values on the horizontal axis correspond to the PDOS 
for the spin majority and minority channels, respectively. The hori- 
zontal dashed lines indicate the Fermi level, (b) Band structure for 
Os/graphene. The solid and dashed lines respectively indicate major- 
ity and minority spin bands, which are widely separated due to strong 
moment formation in the non-relativistic limit, (d) Band structure of 
(Cu-Os)/graphene. In both (b) and (c) the spectrum is always metal- 
lic, demonstrating that the gaps found earlier indeed originate from 
spin-orbit coupling. 
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SUPPLEMENTAL MATERIAL 

Electronic properties without spin-orbit coupling 

To confirm that the band gaps induced by Os and Cu- 
Os stem from spin-orbit coupling, we performed DFT cal- 
culations for these adatoms at 6.25% coverage in the non- 
relativistic limit. For Os/graphene, we find that in the ab- 
sence of spin-orbit coupling the Os spins polarize much more 
strongly compared to the relativistic case yielding a magnetic 
moment M s = 1.52/i^. The spin splitting is clearly visible 
in the partial density of states for the Os bd orbitals shown in 
Fig. [5ja), where positive and negative values on the horizon- 
tal axis correspond respectively to the majority and minority 
spin channels. The band structure for the majority spins (solid 
curves) and minority spins (dashed curves) appears in Fig. 
[5jb). Notice that in sharp contrast to the spin-orbit-coupled 
case, both spin channels are always gapless for any value of 
the chemical potential (at least over the energy range shown). 
Interestingly, (Cu-Os)/graphene remains non-magnetic even 
in the absence of spin-orbit coupling. From the band struc- 
ture in Fig. [5jc) one clearly sees that here, too, the system 
remains gapless in the non-relativistic limit. Thus for both 



Energies of different Cu-Os co-adsorption configurations 

In the main text we remarked that Cu and Os co-adsorbates 
energetically prefer to form vertical dimers over the H-site in 
graphene. Here we summarize the evidence supporting this 
conclusion. We used DFT, without spin-orbit coupling, to ex- 
plore the energetics of the various configurations displayed in 
Fig. [6] (Os, Cu, and C atoms are respectively represented by 
dark cyan, coral, and gray spheres). Note that in all cases 
we fix the Os to the H-site since the binding energy for that 
position greatly exceeds that for the top or bridge sites. (Ad- 
ditionally, Cu-Os dimers bind much more weakly to the top 
and bridge positions.) The vertical Cu-Os dimer in configura- 
tion I exhibits a very large binding energy = 6.17 eV. In 
contrast the binding energy for configurations III and IV are 
substantially lower by 2.50 and 2.09 eV, respectively. Config- 
uration II is unstable and transforms to configuration I after 
relaxation, indicating a strong tendency toward vertical dimer 
formation. As we saw earlier the vertical dimer in (a) is non- 
magnetic (with or without spin-orbit coupling), whereas the 
other configurations are magnetic. 

Results for Ir/graphene 

Finally, we briefly highlight our results for Ir adatoms on 
graphene. We found that Ir slightly prefers the H site over the 
bridge site, the binding energy being 2.17 eV in the former 
configuration and 2.12 eV in the latter. The binding energy 
for Ir at the top site is weakest at 1.94 eV. All results dis- 
cussed henceforth thus correspond to H-site Ir adatoms. With- 
out spin-orbit coupling DFT predicts that Ir/graphene is non- 
magnetic. Figures [TJa) and (b) illustrate the band structure 
and partial density of states for the Ir 5d levels, calculated at 
6.25% coverage. Just as for Os and Cu-Os dimers, here too 
the system remains gapless in the non-relativistic limit. 

Figure [TJd) shows that restoring spin-orbit coupling intro- 
duces a huge band gap A so — 0.21 eV below the Fermi level 
(green dashed line). Moreover, in Fig. [TJc) we see that this 
gap results from hybridization between the Ir d xz / yz orbitals 
and graphene's tt bands, similar to Os/graphene and the tight- 
binding model described in the main text. Spin-orbit interac- 
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FIG. 7: (a) Band structure and (b) PDOS for Ir 5d orbitals obtained 
for Ir/graphene at 6.25% coverage, without SOC. The green hori- 
zontal dashed line denotes the Fermi level. Here the system is non- 
magnetic, and the spectrum is again always gapless similar to Figs. 
[5jb) and (c). In (c) and (d) we display the PDOS and band structure 
calculated with spin-orbit coupling. A large spin-orbit-induced gap 
opens below the Fermi level, though the bands are split slightly due 
to weak moments formed by Ir in this case. As discussed in the main 
text replacing Ir with Cu-Ir dimers returns the Fermi level to within 
the gap and quenches the moment to restore time-reversal symmetry. 



tions also, however, produce small spin and orbital moments 
for Ir given respectively by 0.30 fi b and 0. 13 /jLb \ these are re- 
sponsible for the weak band splittings in Fig.[7Jd). Coverage 
dependence of the gap is illustrated in Fig. [8] As in all other 
cases studied here, the gap remains quite large down to very 
dilute Ir concentrations. Even at 1% coverage where A so 
reduces to 0.08 eV the gap exceeds that induced by In or Tl 
adatoms l23ll many times over. 

Observing a true topological insulator phase in Ir/graphene 
is complicated by the magnetic moments predicted by DFT 
and the fact that the Fermi level resides in the conduction 
band. (One should again keep in mind, however, that DFT 
may overestimate the tendency for moment formation here. 
Furthermore, it may be feasible to gate the system back to the 
insulating regime by conventional means at low coverages.) In 
the main text we described how co-doping with Cu eliminates 
both potential challenges. Here we simply wish to note that 
Cu-Ir dimers carry yet another advantage — they strongly en- 
hance binding to the H-site compared to isolated Ir adatoms, 
similar to the Os case. Indeed, the binding energy for Cu-Ir 
dimers over the H-site is larger by 0.62 eV and 0.79 eV com- 
pared to the bridge and top sites, respectively. By contrast the 
differences for pure Ir are only 0.05 eV and 0.23 eV as we saw 
above. 
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